{smcl}
{com}{sf}{ul off}{txt}{.-}
      name:  {res}<unnamed>
       {txt}log:  {res}/Users/xpmiki/Dropbox/esaiasson/replicationData/e2.smcl
  {txt}log type:  {res}smcl
 {txt}opened on:  {res} 9 Sep 2016, 09:28:45
{txt}
{com}. clear 
{txt}
{com}. set more off
{txt}
{com}. use "symbols teachers.dta"
{txt}
{com}. 
. drop if rep == 1
{txt}(96 observations deleted)

{com}. 
. set seed  358395
{txt}
{com}. sem (d <-  t  wl p)(p <-  t  wl)  
{res}{txt}(70 observations with missing values excluded)

Endogenous variables

{p 0 11 2}Observed:{space 2}{res}d p{p_end}
{txt}
Exogenous variables

{p 0 11 2}Observed:{space 2}{res}t wl{p_end}
{txt}
Fitting target model:

Iteration 0:{space 3}log likelihood = {res:-196.64468}  
Iteration 1:{space 3}log likelihood = {res:-196.64468}  

Structural equation model{col 49}Number of obs{col 67}= {res}       142
{col 1}{txt}Estimation method{col 20}= {res}ml
{col 1}{txt}Log likelihood{col 20}= {res}-196.64468

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}{col 26}      OIM
{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      z{col 46}   P>|z|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}{res}Structural  {col 14}{txt}{c |}
{space 2}{col 3}d <-      {col 14}{c |}
{space 11}p {c |}{col 14}{res}{space 2} .3953561{col 26}{space 2} .0661019{col 37}{space 1}    5.98{col 46}{space 3}0.000{col 54}{space 4} .2657988{col 67}{space 3} .5249134
{txt}{space 11}t {c |}{col 14}{res}{space 2} .0454477{col 26}{space 2} .0392885{col 37}{space 1}    1.16{col 46}{space 3}0.247{col 54}{space 4}-.0315563{col 67}{space 3} .1224517
{txt}{space 10}wl {c |}{col 14}{res}{space 2} .1781897{col 26}{space 2}  .042132{col 37}{space 1}    4.23{col 46}{space 3}0.000{col 54}{space 4} .0956125{col 67}{space 3} .2607668
{txt}{space 7}_cons {c |}{col 14}{res}{space 2} .3753467{col 26}{space 2} .0378617{col 37}{space 1}    9.91{col 46}{space 3}0.000{col 54}{space 4} .3011391{col 67}{space 3} .4495542
{space 2}{txt}{hline 11}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 2}{col 3}p <-      {col 14}{c |}
{space 11}t {c |}{col 14}{res}{space 2} .0741036{col 26}{space 2} .0494886{col 37}{space 1}    1.50{col 46}{space 3}0.134{col 54}{space 4}-.0228924{col 67}{space 3} .1710995
{txt}{space 10}wl {c |}{col 14}{res}{space 2}  .257742{col 26}{space 2} .0489194{col 37}{space 1}    5.27{col 46}{space 3}0.000{col 54}{space 4} .1618617{col 67}{space 3} .3536223
{txt}{space 7}_cons {c |}{col 14}{res}{space 2} .3593783{col 26}{space 2}  .037428{col 37}{space 1}    9.60{col 46}{space 3}0.000{col 54}{space 4} .2860207{col 67}{space 3} .4327359
{txt}{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
     var(e.d){c |}{col 14}{res}{space 2} .0469289{col 26}{space 2} .0055694{col 54}{space 4} .0371896{col 67}{space 3} .0592188
{txt}     var(e.p){c |}{col 14}{res}{space 2} .0756353{col 26}{space 2} .0089763{col 54}{space 4} .0599384{col 67}{space 3} .0954429
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
LR test of model vs. saturated: chi2({res:0})   = {res:     0.00}, Prob > chi2 = {res}     .
{txt}
{com}. estat ic

Akaike's information criterion and Bayesian information criterion

{txt}{hline 13}{c TT}{hline 63}
       Model {c |}        Obs  ll(null)  ll(model)      df         AIC        BIC
{hline 13}{c +}{hline 63}
{ralign 12:.}{col 14}{c |}{res}{col 16}       142{col 27}        .{col 38}-196.6447{col 49}     9{col 58} 411.2894{col 69} 437.8918
{txt}{hline 13}{c BT}{hline 63}
{p 15 21 2}
Note: N=Obs used in calculating BIC; see {helpb bic_note:[R] BIC note}.
{p_end}

{com}. estat teffects 
{res}

{txt}Direct effects
{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}{col 26}      OIM
{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      z{col 46}   P>|z|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}{res}Structural  {col 14}{txt}{c |}
{space 2}{col 3}d <-      {col 14}{c |}
{space 11}p {c |}{col 14}{res}{space 2} .3953561{col 26}{space 2} .0661019{col 37}{space 1}    5.98{col 46}{space 3}0.000{col 54}{space 4} .2657988{col 67}{space 3} .5249134
{txt}{space 11}t {c |}{col 14}{res}{space 2} .0454477{col 26}{space 2} .0392885{col 37}{space 1}    1.16{col 46}{space 3}0.247{col 54}{space 4}-.0315563{col 67}{space 3} .1224517
{txt}{space 10}wl {c |}{col 14}{res}{space 2} .1781897{col 26}{space 2}  .042132{col 37}{space 1}    4.23{col 46}{space 3}0.000{col 54}{space 4} .0956125{col 67}{space 3} .2607668
{space 2}{txt}{hline 11}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 2}{col 3}p <-      {col 14}{c |}
{space 11}t {c |}{col 14}{res}{space 2} .0741036{col 26}{space 2} .0494886{col 37}{space 1}    1.50{col 46}{space 3}0.134{col 54}{space 4}-.0228924{col 67}{space 3} .1710995
{txt}{space 10}wl {c |}{col 14}{res}{space 2}  .257742{col 26}{space 2} .0489194{col 37}{space 1}    5.27{col 46}{space 3}0.000{col 54}{space 4} .1618617{col 67}{space 3} .3536223
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}


Indirect effects
{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}{col 26}      OIM
{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      z{col 46}   P>|z|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}{res}Structural  {col 14}{txt}{c |}
{space 2}{col 3}d <-      {col 14}{c |}
{space 11}p {c |}{col 14}{res}{space 2}        0{col 26}{txt}  (no path)
{space 11}t {c |}{col 14}{res}{space 2} .0292973{col 26}{space 2} .0201695{col 37}{space 1}    1.45{col 46}{space 3}0.146{col 54}{space 4}-.0102342{col 67}{space 3} .0688288
{txt}{space 10}wl {c |}{col 14}{res}{space 2} .1018999{col 26}{space 2} .0257745{col 37}{space 1}    3.95{col 46}{space 3}0.000{col 54}{space 4} .0513827{col 67}{space 3}  .152417
{space 2}{txt}{hline 11}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 2}{col 3}p <-      {col 14}{c |}
{space 11}t {c |}{col 14}{res}{space 2}        0{col 26}{txt}  (no path)
{space 10}wl {c |}{col 14}{res}{space 2}        0{col 26}{txt}  (no path)
{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}


Total effects
{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}{col 26}      OIM
{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      z{col 46}   P>|z|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}{res}Structural  {col 14}{txt}{c |}
{space 2}{col 3}d <-      {col 14}{c |}
{space 11}p {c |}{col 14}{res}{space 2} .3953561{col 26}{space 2} .0661019{col 37}{space 1}    5.98{col 46}{space 3}0.000{col 54}{space 4} .2657988{col 67}{space 3} .5249134
{txt}{space 11}t {c |}{col 14}{res}{space 2}  .074745{col 26}{space 2} .0436166{col 37}{space 1}    1.71{col 46}{space 3}0.087{col 54}{space 4}-.0107419{col 67}{space 3} .1602319
{txt}{space 10}wl {c |}{col 14}{res}{space 2} .2800895{col 26}{space 2} .0431149{col 37}{space 1}    6.50{col 46}{space 3}0.000{col 54}{space 4} .1955859{col 67}{space 3} .3645932
{space 2}{txt}{hline 11}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 2}{col 3}p <-      {col 14}{c |}
{space 11}t {c |}{col 14}{res}{space 2} .0741036{col 26}{space 2} .0494886{col 37}{space 1}    1.50{col 46}{space 3}0.134{col 54}{space 4}-.0228924{col 67}{space 3} .1710995
{txt}{space 10}wl {c |}{col 14}{res}{space 2}  .257742{col 26}{space 2} .0489194{col 37}{space 1}    5.27{col 46}{space 3}0.000{col 54}{space 4} .1618617{col 67}{space 3} .3536223
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}. 
. *** indirect and direct - winning/losing
. cap program drop indireff
{txt}
{com}. program indireff, rclass
{txt}  1{com}. sem (d <-  t  wl p)(p <-  t  wl)  
{txt}  2{com}.   estat teffects
{txt}  3{com}.   mat bi = r(indirect)
{txt}  4{com}.   mat bt = r(total)
{txt}  5{com}.   return scalar indir  = el(bi,1,3)
{txt}  6{com}.   return scalar total  = el(bt,1,3)
{txt}  7{com}. end 
{txt}
{com}. 
. set seed 358395
{txt}
{com}. bootstrap r(indir) r(total), reps(1000): indireff
{txt}(running indireff on estimation sample)

Bootstrap replications ({res}1000{txt})
{hline 4}{c +}{hline 3} 1 {hline 3}{c +}{hline 3} 2 {hline 3}{c +}{hline 3} 3 {hline 3}{c +}{hline 3} 4 {hline 3}{c +}{hline 3} 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000
{res}
{txt}Bootstrap results{col 49}Number of obs{col 67}= {res}       142
{txt}{col 49}Replications{col 67}= {res}     1,000
{p2colset 7 17 21 2}{...}

{txt}{p2col :command:}indireff{p_end}
{p2colset 9 17 21 2}{...}
{p2col :_bs_1:}{res:r(indir)}{p_end}
{p2colset 9 17 21 2}{...}
{p2col :_bs_2:}{res:r(total)}{p_end}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}   Observed{col 26}   Bootstrap{col 54}         Norm{col 67}al-based
{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      z{col 46}   P>|z|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 7}_bs_1 {c |}{col 14}{res}{space 2} .1018999{col 26}{space 2}  .032004{col 37}{space 1}    3.18{col 46}{space 3}0.001{col 54}{space 4} .0391733{col 67}{space 3} .1646265
{txt}{space 7}_bs_2 {c |}{col 14}{res}{space 2} .2800895{col 26}{space 2} .0469054{col 37}{space 1}    5.97{col 46}{space 3}0.000{col 54}{space 4} .1881566{col 67}{space 3} .3720224
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}. estat bootstrap, percentile bc  

{txt}Bootstrap results{col 49}Number of obs{col 67}= {res}       142
{txt}{col 49}Replications{col 67}= {res}      1000
{p2colset 7 17 21 2}{...}

{txt}{p2col :command:}indireff{p_end}
{p2colset 9 17 21 2}{...}
{p2col :_bs_1:}{res:r(indir)}{p_end}
{p2colset 9 17 21 2}{...}
{p2col :_bs_2:}{res:r(total)}{p_end}

{col 1}{text}{hline 13}{c TT}{hline 12}{hline 11}{hline 13}{hline 12}{hline 10}{hline 6}
{col 14}{text}{c |}    Observed{col 38}    Bootstrap
{col 14}{text}{c |}       Coef.{col 27}       Bias{col 38}    Std. Err.{col 51}  [95% Conf. Interval]
{col 1}{text}{hline 13}{c +}{hline 12}{hline 11}{hline 13}{hline 12}{hline 10}{hline 6}
{col 1}{text}       _bs_1{col 14}{c |}{result}{space 2} .10189986{col 27}{space 2} .0011687{col 38}{space 2} .03200395{col 51}{space 2} .0466332{col 63}{space 1} .1758736{col 73}{text}   (P)
{col 14}{text}{c |}{col 51}{result}{space 2} .0509604{col 63}{space 1} .1842398{col 73}{text}  (BC)
{col 1}{text}       _bs_2{col 14}{c |}{result}{space 2} .28008952{col 27}{space 2}-.0001484{col 38}{space 2} .04690539{col 51}{space 2} .1857351{col 63}{space 1}  .370834{col 73}{text}   (P)
{col 14}{text}{c |}{col 51}{result}{space 2}  .183557{col 63}{space 1} .3707279{col 73}{text}  (BC)
{col 1}{text}{hline 13}{c BT}{hline 12}{hline 11}{hline 13}{hline 12}{hline 10}{hline 6}
{col 0}(P){col 8}percentile confidence interval
{col 0}(BC){col 8}bias-corrected confidence interval

{com}. 
. *** indirect and direct - treatment
. cap program drop indireff
{txt}
{com}. program indireff, rclass
{txt}  1{com}. sem (d <- wl  t  p)(p <- wl  t )   
{txt}  2{com}.   estat teffects
{txt}  3{com}.   mat bi = r(indirect)
{txt}  4{com}.   mat bt = r(total)
{txt}  5{com}.   return scalar indirect  = el(bi,1,3)
{txt}  6{com}.   return scalar total  = el(bt,1,3)
{txt}  7{com}. end
{txt}
{com}. 
. set seed 358395
{txt}
{com}. bootstrap r(indirect) r(total), reps(1000): indireff
{txt}(running indireff on estimation sample)

Bootstrap replications ({res}1000{txt})
{hline 4}{c +}{hline 3} 1 {hline 3}{c +}{hline 3} 2 {hline 3}{c +}{hline 3} 3 {hline 3}{c +}{hline 3} 4 {hline 3}{c +}{hline 3} 5 
..................................................    50
..................................................   100
..................................................   150
..................................................   200
..................................................   250
..................................................   300
..................................................   350
..................................................   400
..................................................   450
..................................................   500
..................................................   550
..................................................   600
..................................................   650
..................................................   700
..................................................   750
..................................................   800
..................................................   850
..................................................   900
..................................................   950
..................................................  1000
{res}
{txt}Bootstrap results{col 49}Number of obs{col 67}= {res}       142
{txt}{col 49}Replications{col 67}= {res}     1,000
{p2colset 7 17 21 2}{...}

{txt}{p2col :command:}indireff{p_end}
{p2colset 9 17 21 2}{...}
{p2col :_bs_1:}{res:r(indirect)}{p_end}
{p2colset 9 17 21 2}{...}
{p2col :_bs_2:}{res:r(total)}{p_end}

{res}{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}   Observed{col 26}   Bootstrap{col 54}         Norm{col 67}al-based
{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      z{col 46}   P>|z|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 7}_bs_1 {c |}{col 14}{res}{space 2} .0292973{col 26}{space 2} .0214698{col 37}{space 1}    1.36{col 46}{space 3}0.172{col 54}{space 4}-.0127828{col 67}{space 3} .0713774
{txt}{space 7}_bs_2 {c |}{col 14}{res}{space 2}  .074745{col 26}{space 2} .0439319{col 37}{space 1}    1.70{col 46}{space 3}0.089{col 54}{space 4}-.0113599{col 67}{space 3} .1608499
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}

{com}. estat bootstrap, percentile bc  

{txt}Bootstrap results{col 49}Number of obs{col 67}= {res}       142
{txt}{col 49}Replications{col 67}= {res}      1000
{p2colset 7 17 21 2}{...}

{txt}{p2col :command:}indireff{p_end}
{p2colset 9 17 21 2}{...}
{p2col :_bs_1:}{res:r(indirect)}{p_end}
{p2colset 9 17 21 2}{...}
{p2col :_bs_2:}{res:r(total)}{p_end}

{col 1}{text}{hline 13}{c TT}{hline 12}{hline 11}{hline 13}{hline 12}{hline 10}{hline 6}
{col 14}{text}{c |}    Observed{col 38}    Bootstrap
{col 14}{text}{c |}       Coef.{col 27}       Bias{col 38}    Std. Err.{col 51}  [95% Conf. Interval]
{col 1}{text}{hline 13}{c +}{hline 12}{hline 11}{hline 13}{hline 12}{hline 10}{hline 6}
{col 1}{text}       _bs_1{col 14}{c |}{result}{space 2}  .0292973{col 27}{space 2} -.001192{col 38}{space 2} .02146983{col 51}{space 2}-.0101844{col 63}{space 1} .0760458{col 73}{text}   (P)
{col 14}{text}{c |}{col 51}{result}{space 2}-.0046412{col 63}{space 1} .0852468{col 73}{text}  (BC)
{col 1}{text}       _bs_2{col 14}{c |}{result}{space 2} .07474498{col 27}{space 2}  .000108{col 38}{space 2} .04393186{col 51}{space 2}-.0126474{col 63}{space 1} .1644036{col 73}{text}   (P)
{col 14}{text}{c |}{col 51}{result}{space 2}-.0093307{col 63}{space 1} .1651437{col 73}{text}  (BC)
{col 1}{text}{hline 13}{c BT}{hline 12}{hline 11}{hline 13}{hline 12}{hline 10}{hline 6}
{col 0}(P){col 8}percentile confidence interval
{col 0}(BC){col 8}bias-corrected confidence interval

{com}. 
. *** randomization check
. reg  t  gender interest ideology , robust

{txt}Linear regression                               Number of obs     = {res}       197
                                                {txt}F(3, 193)         =  {res}     0.82
                                                {txt}Prob > F          = {res}    0.4823
                                                {txt}R-squared         = {res}    0.0124
                                                {txt}Root MSE          =    {res} .50186

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}{col 26}    Robust
{col 1}           t{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 6}gender {c |}{col 14}{res}{space 2}  .039834{col 26}{space 2} .0864943{col 37}{space 1}    0.46{col 46}{space 3}0.646{col 54}{space 4}-.1307615{col 67}{space 3} .2104294
{txt}{space 4}interest {c |}{col 14}{res}{space 2}-.0771134{col 26}{space 2} .1413144{col 37}{space 1}   -0.55{col 46}{space 3}0.586{col 54}{space 4}-.3558323{col 67}{space 3} .2016055
{txt}{space 4}ideology {c |}{col 14}{res}{space 2} .1880098{col 26}{space 2} .1381668{col 37}{space 1}    1.36{col 46}{space 3}0.175{col 54}{space 4}-.0845009{col 67}{space 3} .4605205
{txt}{space 7}_cons {c |}{col 14}{res}{space 2}   .47067{col 26}{space 2}  .124931{col 37}{space 1}    3.77{col 46}{space 3}0.000{col 54}{space 4} .2242647{col 67}{space 3} .7170753
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}{txt}
{com}. log close
      {txt}name:  {res}<unnamed>
       {txt}log:  {res}/Users/xpmiki/Dropbox/esaiasson/replicationData/e2.smcl
  {txt}log type:  {res}smcl
 {txt}closed on:  {res} 9 Sep 2016, 09:30:30
{txt}{.-}
{smcl}
{txt}{sf}{ul off}